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We report on the metalorganic chemical vapor deposition (MOCVD) of GaN:Fe and (Ga,Fe)N 
layers on c-sapphire substrates and their thorough characterization via high-resolution x-ray diffrac- 
tion (HRXRD), transmission electron microscopy (TEM), spatially-resolved energy dispersive X-ray 
spectroscopy (EDS), secondary-ion mass spectroscopy (SIMS), photoluminescence (PL), Hall-effect, 
electron-paramagnetic resonance (EPR) , and magnetometry employing a superconducting quantum 
interference device (SQUID). A combination of TEM and EDS reveals the presence of coherent 
nanocrystals presumably Fe x N with the composition and lattice parameter imposed by the host. 
From both TEM and SIMS studies, it is stated that the density of nanocrystals and, thus the 
Fe concentration increases towards the surface. According to Hall effect measurements, electrons 
from residual donors are trapped by mid-gap Fe acceptor states in the limit of low iron content 
x < 0.4%, indicating that the concentration of Fe 2+ ions increases at the expense of Fe ions in the 
3+ charge state. This effect is witnessed by photoluminescence (PL) measurements as changes in 
the intensity of the Fe 3+ -related intra-ionic transition, which can be controlled by co-doping with Si 
donors and Mg acceptors. In this regime, EPR of Fe 3+ ions and Curie-like magnetic susceptibility 
are observed. As a result of the spin-orbit interaction, Fe 2+ does not produce any EPR response. 
However, the presence of Fe ions in the 2+ charge state may account for a temperature-independent 
Van Vleck-type paramagnetic signal that we observe by SQUID magnetometry. Surprisingly, at 
higher Fe concentrations, the electron density is found to increase substantially with the Fe content. 
The co-existence of electrons in the conduction band and Fe in the 3+ charge state is linked to the 
gradient in the Fe concentration. In layers with iron content x > 0.4% the presence of ferromagnetic 
signatures, such as magnetization hysteresis and spontaneous magnetization, have been detected. A 
set of precautions has been undertaken in order to rule out possible sources of spurious ferromag- 
netic contributions. Under these conditions, a ferromagnetic- like response is shown to arise from 
the (Ga,Fe)N epilayers, it increases with the iron concentration, it persists up to room temperature, 
and it is anisotropic - i.e., the saturation value of the magnetization is higher for in- plane magnetic 
field. We link the presence of ferromagnetic signatures to the formation of Fe-rich nanocrystals, 
as evidenced by TEM and EDS studies. This interpretation is supported by magnetization mea- 
surements after cooling in- and without an external magnetic field, pointing to superparamagnetic 
properties of the system. It is argued that the high temperature ferromagnetic response due to 
spinodal decomposition into regions with small and large concentration of the magnetic component 
is a generic property of diluted magnetic semiconductors and diluted magnetic oxides showing high 
apparent Curie temperature. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.50.Tt, 81.05.Ea 



I. INTRODUCTION 

In recent years, it has become more and more clear that 
wide band-gap semiconductors and oxides doped with 
transition metals (TM) constitute a new class of mate- 
rial systems exhibiting magnetic properties whose origin 
and methods of control are still not understoodi^^^ 
To this category belongs certainly (Ga,Fe)N, which is 
the subject of the present work. While extensive stud- 
ies have been conducted on (Ga,Mn)N^^^i as promising 
workbench for future applications in spintronics^^ only 
little is known about the (Ga,Fe)N materials system. 



The magnetic properties of nominally undoped and p- 
doped GaN implanted with Fe ions were reported, show- 
ing a ferromagnetic response evidenced by magnetiza- 
tion hysteresis loops persisting, depending on the pro- 
vided dose,— up to room temperature (RT) . 11 ' 12 Fur- 
thermore, the emission channeling technique applied to 
Fe-implanted (Fe dose up to 10 19 cm" 3 ) GaN samples 
confirmed the presence of a high percentage of TM ions 
(up to 80%) occupying substitutional Ga sites of the 
host crystalfi 3 - Moreover, GaN films doped with Fe, with 
concentrations up to 3 x 10 19 cm -3 were grown by 
molecular beam epitaxy (MBE) at substrate tempera- 
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tures Tg ranging from 20 to 380° C directly on sapphire 
(0001) and ferromagnetic behavior with a Curie temper- 
ature Tq = 100 K has been observed only in the sam- 
ples grown at 400°C,H Films of GaN:Fe (Fe density 
up to 6 x 10 21 cm~ 3 ) fabricated by means of MBE 
at Ts = 500 - 800° C showed a superparamagnetic 
behavio r 15 ' 16 assigned to Ga-Fe and/or Fe-N inclusions. 
Extended x-ray absorption fine structure (EXAFS) anal- 
ysis suggests that the decrease of Tq leads to a struc- 
tural transition from the wurtzite (wz) to the zinc-blende 
(zb) reconstruction, and this transition may be related to 
the origin of ferromagnetism in Fe-doped GaN films. 17 
The metalorganic chemical vapor deposition (MOCVD) 
of GaN:Fe has been previously reported, with a focus 
on the actual Fe content in the layers and its effect 
onto the carrier concentration. 18 More recent prelimi- 
nary works by others^ and u o 20 ' 21 indicate that MOCVD 
grown (Ga,Fe)N shows ferromagnetic-like characteristics 
to above room temperature. Moreover, diluted magnetic 
semiconductors, and (Ga,Fe)N in particular, have be- 
come model systems to test various implementations of 
the density functional theory to disordered strongly cor- 
related systems ^^24^5 

Despite the mentioned previous reports, a thorough 
characterization of this novel materials system, with spe- 
cial attention to the mechanisms of the magnetic re- 
sponse, is still needed. In order to predict the feasibility 
of e. g. carrier-mediated spin coupling, it is crucial to 
know the electronic and spin structure of the TM cen- 
ters, which depends on the actual charge state and is 
strongly influenced by the host crystal. Concerning the 
position of the deep acceptor- like level Fe 3+ /Fe 2+ in the 
GaN gap, it was estimated from photoluminescence (PL) 
excitation spectroscopy studies to be 2.5 - 2.6 eV above 
the GaN valence band edge E v . 2e ^ 27 However, another 
group suggests that this level may lie much closer to the 
conduction band at (E v + 2.863) ± 0.005 eV^ 

The present work is devoted to a comprehensive study 
of the GaN:Fe (below the solubility limit of Fe into GaN) 
and (Ga,Fe)N material systems, beginning with a careful 
on-line control of the growth process and proceeding with 
a possibly thorough investigation of the structural, elec- 
trical, optical, and magnetic properties in order to shed 
new light into the mechanisms responsible for the para- 
magnetic and high-temperature ferromagnetic response 
of these novel systems. Particular attention is paid to 
avoid contamination of the superconducting quantum 
interference device (SQUID) signal by spurious effects 
and, thus, to reliably determine the magnetic proper- 
ties of (Ga,Fe)N. By combining various characterization 
techniques, we establish experimentally the relation be- 
tween structural, electronic, and magnetic properties of 
the studied material. We show how doping over the solu- 
bility limit affects the sample morphology and magnetic 
ion distribution, which in turn affect in a dramatic way 
both electronic and magnetic properties of the system. 
The interpretation of our results in terms of spinodal de- 
composition and non-uniform distribution of magnetic 



nanocrystals along the growth direction can be applied 
to explain ferromagnetic-like properties persisting up to 
high temperatures in a broad class of diluted magnetic 
semiconductors and diluted magnetic oxides. 

Our paper is organized as follows. In the next sec- 
tion we describe the experimental methods employed 
to characterize in-situ and ex-situ the structural, opti- 
cal, transport and magnetic properties of the MOCVD 
structures under investigation. Section III contains a 
summary of sample parameters, including Fe concentra- 
tions, as determined by various experimental methods. 
In Sec. IV we present experimental results starting from 
findings obtained from transmission electron microscopy 
(TEM) and spatially-resolved energy dispersive X-ray 
spectroscopy (EDS). These measurements reveal a highly 
non-uniform distribution of the magnetic ions at Fe con- 
centrations surpassing the solubility limit. We then dis- 
cuss how the presence of magnetic nanocrystals and the 
non-uniformity in their distribution modify the electrical, 
optical, magnetic resonance, and magnetic characteris- 
tics of the system. Conclusions and outlook stemming 
from our work are summarized in Sec. V. Finally, in Ap- 
pendix, we discuss in detail the procedure employed to 
determine reliably the magnetic response of (Ga,Fe)N by 
SQUID magnetometry. 



II. EXPERIMENTAL 



A. Growth procedure 

The studied epilayers have been fabricated in an 
AIXTRON 200RF horizontal-tube MOCVD reactor. 
All structures have been deposited on c-plane sap- 
phire substrates according to a well established growth 
procedure 2 ^ involving TMGa, NH 3 , Cp 2 Mg, SiH 4 , and 
Cp2Fe as precursors for respectively Ga, N, Mg, Si, and 
Fe, with H2 as a carrier gas. Upon the nitridation of 
the substrate in the reactor, the deposition of a low- 
temperature (540 °C) GaN nucleation layer (LT-NL), its 
annealing under NH3 and the growth of a 1 /im thick 
device- quality GaN layer at 1050°C, as a first step, sev- 
eral series of 500 nm thick GaN:Fe at different substrate 
temperatures (ranging from 750 to 1050°C) and different 
Cp2Fe flow rates [between 50 and 400 standard cubic cen- 
timeters per minute (seem)] have been fabricated. The 
nominal Fe content in subsequently grown samples has 
been alternatively switched from low to high and vice- 
versa, in order to minimize long term memory effects 
due to the presence of residual Fe in the reactor^ The 
samples have been continuously rotated during growth to 
promote the deposition homogeneity. 
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FIG. 1: In-situ ellipsometry: ellipsometric angle ^ at 2 eV as 
a function of time for a standard deposition sequence. Par- 
ticular regions relative to the growth process and marked (a) 
to (g) are discussed in the main text. 



2. In-situ X-ray diffraction 

Furthermore, the growth process has been routinely 
monitored in-situ by x-ray diffraction. In addition to 
the optical windows necessary for the on-line ellipsome- 
try characterization, the MOCVD reactor has been im- 
plemented with two Be windows transparent to the radi- 
ation energy provided by a PANalytical Cu x-ray tube. 
A Johansson monochromator focuses the beam on the 
growing crystal in a geometry suitable for the observa- 
tion of the (1124) reflection of wz-GaN. The diffracted 
beam is collected by a multichannel detector PANalyti- 
cal X'Celerator with an angular resolution of 0.013° per 
pixel in this geometry. This diffraction set-up allows us 
to perform kinetic measurements yielding information on 
the layer composition, thickness and crystalline quality 
even on rotating samples, thanks to a recently devel- 
oped wobbling compensation algorithm i 34 i 35 Moreover, 
real time reciprocal space mapping can be carried out in 
a wide temperature range (20 - 1050°C), though limited 
to a static sample configuration, in order to eliminate 
wobbling effects^ 



B. In-situ monitoring 



1. Spectroscopic ellipsometry and laser inter ferometry 



The employed MOCVD reactor is equipped with an 
in-situ Isa Jobin Yvon ellipsometer which allows both 
spectroscopic (optical response - in terms of the ellipso- 
metric angles 'F and A - as a function of the photon en- 
ergy) and kinetic (optical response as a function of time) 
on-line measurements over the energy range between 1.5 
and 5.5 eV for process optimization.— Figure [1] shows the 
ellipsometric angle $ at 2 eV as a function of time for 
the standard growth process employed to fabricate the 
samples under investigation: (a) heating of the sapphire 
substrate from RT to 1200°C; (b) desorption under H2 
flow to stabilize the surface; (c) cooling and nitridation of 
the substrate; (d) deposition of the LT-NL; (e) anneal- 
ing of the NL; (f) growth of the GaN buffer layer; (g) 
growth of the Fe doped GaN layer. Spectroscopic ellip- 
sometry measurements carried out in real time allow the 
on-line control over the growth-rate and the composition 
of the growing layer and provide additional qualitative 
information on the surface roughness. Moreover, simul- 
taneously to ellipsometry, standard on-line dynamic op- 
tical reflectivity 3 ^, at 2 eV has been routinely performed, 
yielding analogous information on the growing layer as 
ellipsometry, but being more sensitive to wobbling (un- 
avoidable during the MOCVD process) and temperature 
effects^ 3 - 
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FIG. 2: In-situ XRD kinetic measurement: FWHM of the 
(1124) GaN reflection vs. deposition time for two samples 
grown with different Cp2Fe flow rates. The evolution of the 
FWHM during the growth of the GaN buffer is reproducibly 
independent of the sample, whereas the broadening during 
the deposition of the Fe-containing layer is a function of the 
magnetic ions content. 

Kinetic in-situ XRD measurements were routinely car- 
ried out during growth. A full-width-at-the-half maxi- 
mum (FWHM) analysis of the GaN peak in the in-situ 
spectra as a function of the growth time, yields infor- 
mation on the crystal quality of the layers. In Fig. [5] 
a clear decrease of the FWHM during the deposition of 
the GaN buffer due to the growing thickness of high crys- 
talline quality GaN is detectable and was proven to be 
quantitatively reproducible on a full series of samples. 34 
The vertical line demarcates the time corresponding to 
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the opening of the Cp2Fe source and indicates the onset 
of broadening of the GaN peak due to the incorporation 
of Fe ions. A clear dependence of the FWHM on the 
layer thickness and on the iron precursor flow rate is evi- 
dent and comparable to the values obtained from ex-situ 
post-growth XRD analysis reported in Table I. 
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FIG. 3: In-situ XRD space maps in the (1124) reflection of 
GaN and broadening in the inverse of the lattice parame- 
ter along the growth direction (left panels) for two values of 
Cp2Fe flow rates, 50 and 350 seem (upper and lower panel, 
respectively) . 



C. Ex-situ characterization methods 

1. High-resolution X-ray diffraction 

High-resolution X-ray diffraction (HRXRD) rocking 
curves of the GaN (0002) reflex, give for all the sam- 
ples, including the (Ga,Fe)N structures, a FWHM of 
200 - 320 arcsec, comparable with state-of-the-art device- 
quality gallium nitride material^ The broadening of the 
GaN(0002) reflex is found to depend both on the iron 
doping level and on the growth temperature. The mini- 
mum FWHM is obtained for the lowest Cp2Fe flow rate, 
50 seem, and the highest growth temperature, 950°C. 
Despite the importance of HRXRD for the routine struc- 
tural characterization of the layers, the inadequacy of 
standard X-ray diffraction methods for nanoscale inves- 
tigations (detection of the presence of non-uniform mag- 
netic ion distribution) is already widely accepted^ We 
also note that, according to atomic force microscopy 
( AFM) , the surface roughness increases with the Fe con- 
centration, an effect that may affect some quantitative 
considerations. The change in sample morphology indi- 
cates the onset of a three-dimensional (3D) granular-like 
growth mode characterized by grain height differences 
exceeding the layer nominal thickness at the highest Fe 
concentration. 



A more detailed information on the structural proper- 
ties of the layers, can be obtained by acquiring in-situ 
reciprocal space maps. The diffraction geometry for the 
observation of the (1124) reflection has been chosen in 
such a way that the sensitivity to the lattice parameter 
Q z perpendicular to the (0001) surface could be opti- 
mized. After the deposition of the GaN buffer the sub- 
strate rotation was interrupted at the proper azimuth 
for the acquisition of reciprocal space maps. Figure [3] 
presents the space map for the GaN buffer layer (GaN 
@ 850°C) qualitatively reproducible for all investigated 
samples. The procedure was repeated just after the depo- 
sition of the Fe-doped layer (GaN:Fe @ 850°C) and after 
cooling the sample to room temperature (GaN:Fe @ RT) 
upon completion of the growth process. The reciprocal 
maps in Fig. [3] evidence the variation of the lattice pa- 
rameter along the growth direction, AQ Z , obtained from 
the following relation^ 



2. Secondary-ion mass spectroscopy 

The total Fe concentration in the layers has been deter- 
mined via secondary-ion mass spectroscopy (SIMS) and 
calibrated with an undoped GaN sample implanted with 
10 16 cm~ 2 56 Fe at 270 keV as reference. The total Fe 
concentration in the layers is found to increase with the 
nominal Cp2Fe flow rate. We also detect an increase in 
the oxygen concentration which reaches 10 20 cm -3 at the 
highest Cp2Fe flow rates, an effect that may result from 
a larger surface area due to the above-mentioned onset 
of granular growth. We underline that the granular mor- 
phology may render the quantitative determination of Fe 
and O contents and profiles not reliable in the case of the 
highest Fe concentration. 



AQ Z = — • cos(e GaN ) • Ae; (1) 

where A is the wavelength of the incident radiation, ecaN 
gives the diffraction angle for relaxed GaN and Ae is the 
measured relative peak position. As confirmed by in- 
situ kinetic measurements and post growth analysis, the 
broadening is a function of the iron flow rate. 



3. Electron-paramagnetic resonance 

Electron paramagnetic resonance (EPR) data have 
been acquired with a BRUKER ELEXSYS E-580 spec- 
trometer, operating at a frequency of 9.48 GHz (X-band). 
The spectrometer is equipped with a continuous- flow Ox- 
ford cryostat, enabling measurements in the temperature 
range 2 - 300 K. 
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4- Photoluminescence 

The photoluminescence (PL) measurements in ultra- 
violet (UV) regime have been carried out in a set-up 
consisting of a Helium Cadmium (HeCd) laser with a 
wavelength of 325 nm (3.82 eV) and an excitation power 
of 40 mW, a Jobyn Yvon monochromator with a focal 
length of 550 mm and a Jobyn Yvon 1024 x 512 pixel liq- 
uid nitrogen-cooled Charge Coupled Device (CCD) cam- 
era for detection. For the measurements in the infrared 
(IR) region the 351.1 nm (3.53 eV) line of a Spectra 
Physics argon-laser, a Spex Spectrometer with 850 mm 
focal length and a liquid nitrogen cooled Ge-detector were 
employed. In both set-ups, the samples were placed in a 
He-cooled flow cryostat, which allows to cool the samples 
to 5 K. For adjusting the temperature a silicon diode and 
a Lakeshore 331 temperature controller were employed. 



5. Magnetic characterization 

The magnetic properties have been investigated in a 
home made 1 T SQUID magnetometer working in the 
temperature range 5 - 330 K. The samples for magnetic 
investigations are typically cut into 9 x 3 or 5 x 5 mm 2 
pieces, for in-plane and both in- and out-of-plane com- 
parative measurements, respectively. We study both the 
temperature dependence of the magnetization at a con- 
stant field and the sample response to the variation of the 
external field at a constant temperature. In this case, we 
find it sufficient to sweep the field in one direction only - 
a half hysteresis loop - and to obtain the second half by 
means of numerical reflection of the data. Selected mea- 
surements of the full hysteresis loop have been carried 
out in order to ensure the reliability of this experimental 
procedure. Also, all the magnetization data presented in 
the paper are already corrected for the diamagnetic signal 
from the substrate, reflecting, in this way, the properties 
of the investigated layers only. A detailed discussion of 
the substraction procedure and related experimental dif- 
ficulties specific to magnetic measurements on thin films 
of diluted magnetic semiconductors and the precautions 
undertaken in the present work are outlined in the Ap- 
pendix. 



6. Transmission electron microscopy 

Transmission electron microscopy studies have been 
carried out on cross-sectional samples prepared by stan- 
dard mechanical polishing followed by Ar + ion milling 
at 4 kV for about 1 h. Conventional diffraction contrast 
images in bright-field imaging mode and high-resolution 
phase contrast pictures were obtained from a JEOL 2011 
Fast TEM microscope operating at 200 kV and capable of 
an ultimate point-to-point resolution of 0.19 nm and al- 
lowing to image lattice fringes with a 0.14 nm resolution. 



The energy dispersive x-ray analysis has been performed 
via an Oxford Inca EDS equipped with a silicon detector. 



III. SAMPLE PARAMETERS 

In Table I we summarize the data related to the 
GaN:Fe, (Ga,Fe)N and reference GaN samples consid- 
ered in this work. As already mentioned, all samples 
have been grown on c-sapphire substrates. The Fe-doped 
structures consist of a nominally 500 nm thick GaN:Fe or 
(Ga,Fe)N layer deposited onto a 1 thick GaN buffer, 
whereas the reference samples are mere 1 /jm thick GaN 
buffers (nominally undoped) on sapphire. In the table, 
the following values are listed: the Cp2Fe flow rate em- 
ployed to grow the TM-doped layers, the FWHM of the 
(0002) reflex from GaN determined by ex-situ HRXRD, 
the Fe 3+ concentration according to SQUID data and 
evaluated from the Curie constant (in the case of low 
Fe concentration, where the low-temperature paramag- 
netic contribution is not covered by the ferromagnetic 
signal) and from the saturation value of the ferromag- 
netic magnetization, respectively. The sample numbers 
in the figures throughout the paper refer to Table I. 

IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Transmission electron microscopy 

According to our TEM studies, in the limit of low 
Fe-doping, i.e., for Cp2Fe flow rates in the range 50 - 
150 seem the structure of the GaN matrix is not de- 
tectably affected. As indicated in Fig.[3J the lattice spac- 
ing between the (0002) planes is 0.250 ± 0.010 nm, per- 
fectly matching that of undoped GaN. In addition, the 
EDS spectra do not point to Fe aggregation (Fig. [5]), indi- 
cating that the Fe ions are randomly distributed over the 
lattice, presumably in substitutional sites of Ga. Above 
175 seem, however, TEM studies reveal the presence of 
nanocrystals embedded in the host GaN. 

According to the TEM micrographs presented in 
Figs.[Ha) and (b), the precipitates are located both close 
to the surface and inside the epilayer in dislocation-free 
regions, suggesting that precipitation is not related to the 
presence of dislocations. On the other hand, Figs. [6fc) 
and (d) prove that some nanocrystals are embedded into 
the upper surface layers, while some others are decorat- 
ing the screw- or mixed- type threading dislocations in 
a volume proximal to the surface. Furthermore, EDS 
spectra demonstrate that the concentration of Fe ions 
is largely enhanced within the nanocrystals. As shown 
in Fig. [7J the Fe-K a i peak at 6.404 KeV can be seen 
clearly in the spectrum taken in the region of a precipi- 
tate, while in the adjacent areas, the presence of Fe is not 
detected. It should be mentioned that the precipitates 
can be found only in the upper region of the layer within 
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TABLE I: Data related to the investigated GaN, GaFe:N and (Ga,Fe)N samples. The following values are listed: the Cp2Fe 
flow rate employed to grow the Fe-doped layers, the FWHM of the (0002) reflex from GaN determined by ex-situ HRXRD, 
the Fe concentration according to the Curie constant from SQUID data (in the case of low Fe concentration, where the 
Curie component is not covered by the ferromagnetic signal; S =5/2 and g = 2 is assumed) and from the saturation value of 
ferromagnetic magnetization Ms, assuming S = g — 1. 
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FIG. 4: HRTEM lattice fringe image of the (0002) planes 
demonstrating that low Fe-doping (sample #465) does not 
cause detectable strain, the lattice parameter (dooo2 = 0.25 ± 
0.01 nm) being the same as that of undoped GaN. 



300 nm from the surface, and thus well above the inter- 
face between the nominally undoped GaN buffer and the 
Fe-doped layer. The segregation of TM-rich nanocrystals 
towards the surface of the samples, seems to represent a 
general phenomenon in highly TM-doped MOCVD III-V 
compounds: a similar effect has been observed for MnAs 
nanocrystals in MOCVD grown (Ga,Mn)As. 38 

We estimate that the nanocrystals occupy less than 
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KeV 

FIG. 5: EDS spectrum showing the presence of Fe in the 
measured regions of the GaN layers doped with a low concen- 
tration of Fe ions (sample #465) . 



1% of the film volume and their diameter varies from 5 
to 50 nm. Some precipitates are elongated in the growth 
direction, as shown in Fig. EJa), such a behavior being 
expected within the spinodal decomposition scenario-^. 
This non-spherical shape may lead to an additional mag- 
netic anisotropy that can elevate the blocking temper- 
ature. In Fig. Ob) the selected area diffraction pattern 
(SADP) along the (1010) zone axis taken around the pre- 
cipitate in Fig. [Ha) is shown. It reveals the presence 
of an additional hexagonal phase with almost the same 
orientation as the GaN matrix. In the corresponding 
schematic indexed pattern (Fig. 0(c)) closed circles cor- 
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(c) 




FIG. 6: Bright-field transmission electron micrographs of (Ga, 
Fe)N demonstrating the presence of precipitates as indicated 
by the white triangles (a) close to the surface and (b) in the 
mid of the epilayer. Some precipitates (c) are embedded into 
the surface and others, (c) and (d) sink to the dislocations. 



£-Fe3N is formed in the considered precipitates. How- 
ever, it is also expected that the GaN host can stabilize 
Fe x N or (Ga,Fe)a;N nanocrystals in a form not existing 
in the case of nanostructured Fe-N material systems. Ac- 
tually, such a situation has recently been revealed in the 
case of MBE grown (Ge,Mn), where spinodal decompo- 
sition into Ge and novel Mn-rich (Ge,Mn) nanocrystals 
has been found4^ 
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FIG. 8: (a) Elongated precipitate observed in an HRTEM 
image of Moire fringes contrast, (b) SADP pattern acquired 
in the region around the precipitate along the (1010) zone 
axis; (c) the corresponding schematic graph for indexing of 
the diffraction spots. 
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FIG. 7: EDS spectra for the sample #466 (left panel) taken 
around the precipitate (right panel) showing that the con- 
centration of Fe is significantly enhanced in the region of the 
precipitate (recorded as spectrum 3) as compared with the 
surrounding matrix (recorded as spectrum 4) . 
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respond to GaN, open circles to the additional phase and 
crosses to the double diffraction between them. The lat- 
tice parameters of this new phase are a = 0.268 ± 0.010 
nm, c = 0.436 ± 0.010 ran. By comparing the lattice 
constants to those of Fe and iron nitride nanostructures 
and thin films listed in Table ITTl we tend to believe that 



FIG. 9: (a) HRTEM image of a precipitate residing on the 
surface of an highly Fe-doped sample (#466); (b) the corre- 
sponding SADP pattern along the (1010) zone axis. 

Our HRTEM studies allow to detect also precipitates 
residing at the surface of the highly Fe-doped samples, as 
shown in Fig.^a). The SADP pattern in Fig.^b) sug- 
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TABLE II: Structural and magnetic parameters of some iron 
and iron-nitride phases in nanostructures and thin films. To 
be compared with our experimental values. 

lattice parameter 



structure a(nm) c(nm) /ib 



FeN 

7'-Fe 4 N 
7 "-FeN .9i 
e - Fe 3 N 

7 "'-FeNo. 5 -o.7 
? - Fe 2 N 
7-Fe 
a-Fe 



ZB 
RS 



0.430 a 
0.40 6 



AF or 0" 



RS 0.382 6 2.21 a 

ZB 0.433" 1.7 or a 

WZ 0.270 c 0.433 c 1.9 d 

RS 0.450 a ferro" 

HCP 0.443" 0.484 a para/1. 5 h 
ZB 0.361-0.370 6 0.3-2. e 

BCC 0.286 / 2.2 9 



"Reference [41j and references therein 

'Reference 43 

c Reference 4^ 

^Reference 43 

e Reference 4a 

^Reference 4a 

9 Reference 43 

^Reference 4a 



B. Electrical properties 

It has been found previously that Fe doping reduces 
the electron concentration in GaN and leads to a semi- 
insulating material i 49 i 50 i 51 In order to find out how the 
electronic properties evolve when the Fe concentration 
surpasses the solubility limit, we have undertaken electri- 
cal and optical measurements. Together with EPR spec- 
troscopy, such studies allow us to assess the charge and 
spin state of Fe ions controlling the magnetic properties 
of the material system under investigation. Hall effect 
and resistance measurements have been carried out in 
the Van der Pauw geometry on 5 x 5 mm 2 square cuts of 
GaN, GaN:Fe, and (Ga,Fe)N wafers, whose characteris- 
tics are summarized in Table I. Soldered indium contacts 
are employed and the ohmic characteristics of all contacts 
are carefully tested. The Hall resistance has been mea- 
sured in the magnetic field B — ±500 mT using currents 
in the range 10 to 200 /iA. The employed current values 
are high enough to achieve high Hall voltages for better 
accuracy, but still sufficiently low to avoid Joule heating 
of the sample. 

In the studied samples, both the (Ga,Fe)N film and 
the GaN buffer layer are conducting, so that a procedure 
suitable for the evaluation of electrical characteristics in 
the presence of a parallel conducting channel has to be 
employed. Accordingly, in order to obtain the conduc- 
tance tensor components G,j of the (Ga,Fe)N films, the 
buffer contribution has to be subtracted from the 



total conductance G; 7 -. We evaluate and G) 1 ,' from 
the measured values of the square resistance R and Hall 
resistance i?Haii for the reference GaN sample and for the 
samples containing the (Ga,Fe)N films according to: 



G%=t b RV/[t r (RV +< u )] 



(2) 



gests that the Moire fringes in this precipitate (Fig.[9^a)) 
are formed by the interference of the GaN (0002) planes 
of the matrix crystal with a set of lattice planes parallel to 
(0002) with a d-spacing of 0.204 ± 0.010 nm, which is al- 
most the same as that of the (111) planes of 7-Fe or of the 
(110) planes of a-Fe, as reported in Table HT1 Therefore, 
these precipitates segregating to the surface could con- 
sist of pure iron, this assumption being consistent with 
the evaporation of nitrogen from the surface during the 
MOCVD process and with the consequent hindering of 
the Fe-nitride compounds formation. 

One may expect that the highly non-uniform Fe distri- 
bution over the layer volume shown by TEM investiga- 
tions will have a dramatic influence on the electronic and 
magnetic properties of (Ga,Fe)N. Experimental results 
presented in the next sub-sections confirm this conjec- 
ture. 



G ( yl = t b R£l u /[t r (R^ 2 + 4 r i u 2 )]; (3) 
G XX = R/[{R 2 + Rl &n )]; (4) 

G yx = i?Hall/[(-R 2 + #Hall)]; ( 5 ) 

where t, tf,, and t r are the thickness of the (Ga,Fe)N 
film, of the buffer layer, and of the reference sample, 
respectively. Therefore, the (Ga,Fe)N conductance is 
= Gij — Gy and from that we obtain the resistivity 
p and the Hall electron concentration n of the (Ga,Fe)N 
films in the standard way: 

p = t/GiQ(B = 0); (6) 

n = B(GiQ 2 + Gif)/(etGgy). (7) 

We must point out that within this model we assume 
the carrier concentration of the GaN buffer layer to be 



9 



the same in all the measured samples and that there is 
no diffusion of Fe into the buffer. Fluctuations in the 
buffer properties will affect the data in a significant way 
when the conductance of the (Ga,Fe)N layer becomes 
small. All the studied samples are found to be n-type 
with Rh < 0, where Rh = RHair t/B is the Hall constant. 
As shown in Fig. [TOJ in the region of low Cp2Fe flow rate 
we find the magnitude and temperature dependencies of 
the Hall concentration typical for MOCVD GaN epitax- 
ial films, in which electrons freeze-out at residual donors 
at low temperatures. The obtained activation energy of 
the carriers in GaN, namely E a = 31 meV, supports 
this conclusion. According to our findings summarized 
in Figs. [TOl and fTTl the samples grown at the low Cp2Fc 
flow rate of 150 seem show reduced Hall concentration as 
compared to GaN, which may result from trapping of the 
electrons by Fe ions, Fe 3+ — > Fe 2+ , an effect consistent 
with the position of the Fe 3+ /Fe 2+ acceptor level in the 
band gap of GaN, as we will discuss in more detail in the 
next subsection. 
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FIG. 10: Electron concentration n from Hall effect measure- 
ments as a function of the inverse temperature for a GaN layer 
and a series of GaN:Fe - (Ga,Fe)N samples with differing Fe 
content controlled by the Cp2Fe flow rate. 



Interestingly, when the Cp2Fe flow rate exceeds 
150 seem, which according to SIMS data corresponds to a 
Fe content above 10 20 cm~ 3 , we observe a strong increase 
of the Hall concentration, and a change in the character 
of its temperature dependence, a behavior clearly seen in 
Figs. [TOl and[TTl For the highest flow rate of 350 seem, the 
Hall concentration reaches 10 18 cm -3 , a value more than 
two orders of magnitude higher than that for the sample 
grown with the Cp2Fe flow rate of 150 seem. This effect 
might suggest the appearance of an additional conduc- 
tance channel associated with the hopping of electrons 
between Fe ions at large magnetic ions densities. Ac- 
tually, hopping conductance involving Fe impurities was 
discussed for InP, 5 - where, however, the measured re- 
sistance was found to be much larger than in our case. 
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FIG. 11: Electron concentration n from Hall effect measure- 
ments as a function of the Cp2Fe flow rate for two GaN ref- 
erence layers and for two series of (GaN:Fe - Ga,Fe)N sam- 
ples (triangles and squares, respectively) at room tempera- 
ture. Reduced electron concentrations in the low flow rate 
region points to a semi-insulating character of Fe-doped lay- 
ers in this regime. 



Indeed, the relatively large values for the Hall mobility 
fi = l/(enp) shown in Fig. 1121 make the hopping scenario 
improbable. We also note that, as discussed in subsec- 
tion |IV3 Fe-rich nanocrystals are formed at higher flow 
rates, so that we deal with carrier transport in a highly 
inhomogenous medium. Under such conditions, the rela- 
tion between the Hall and carrier densities becomes intri- 
cate. It is more probable that (at least a part) of the rise 
in the Hall concentration is caused by an increase in the 
oxygen donor density for large Cp2Fe flow rates. SIMS 
analysis shows that the oxygen concentration increases 
by over two orders of magnitude upon changing the flow 
rate from 50 to 350 seem and at the highest flow rate it 
becomes comparable to the iron concentration. Finally, 
it is also possible that some part of the Fe ions occu- 
pies interstitial positions in which Fe 2+ acts as a double 
donor. Optical and ESR data presented in the next sub- 
sections show that the charge state of Fe ions transforms 
from +3 to +2 with increasing Fe concentration. This 
effect supports the scenario which assigns the enhance- 
ment in the carrier concentration to the increase of the 
donor concentration with the Fe content, though the ori- 
gin of the relevant impurity or defect remains uncertain. 
We will return to this open question when discussing the 
results of the EPR measurements. 



C. Optical properties 

Photoluminescence measurements in the infrared, visi- 
ble, and ultra-violet regime have been carried out on the 
whole series of samples with increasing Fe concentration. 
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FIG. 12: Electron Hall mobility (i as a function of tempera- 
ture for a GaN layer and a series of (Ga,Fe)N samples with 
differing Fe content controlled by the Cp2Fe flow rate. 



FIG. 13: Photoluminescence spectra at 10 K: comparison be- 
tween samples with low (#409) and high (#466) Fe concen- 
tration. 



Figure [13] presents PL spectra acquired in the energy 
range 2.0 - 3.6 eV at 10 K for two GaN:Fe samples grown 
with iron precursor flow rates of 50 and 350 seem. In both 
samples a weak defect-related yellow luminescence band 
(YL) centered around 2.2 eV is seen, as well as the exci- 
tonic near-band-edge emission. The intensity of the latter 
is found to increase by four orders of magnitude between 
the samples with the lowest and the highest Fe content, 
as shown in Figure [141 The peak position varies between 
3.4811 and 3.4833 eV at 10 K, depending on the layer 
thickness. The determined position and thermal behav- 
ior is identical to that of the neutral-donor-bound exciton 
D°X J 4, related to the upper valence subband of wz-GaN 
and observed in nominally undoped GaNi^ 3 - However, the 
peak is noticeably wider, with a FWHM of 4.65 meV for 
the lowest iron concentration and broadens with increas- 
ing doping level. 

The lower intensity emission at 3.4609 eV visible in 
Fig. [15] is assigned to a two electron satellite D X^„ = 2, 
where the recombination of the D°Xa leaves the donor in 
its excited n = 2 statei 54 ' 55 From the separation between 
the principal D°X J 4 and the D°X J 4 i „ = 2 maximum, the 
donor-ionization energy can be obtained^ In our case it 
results in ss 30 meV, pointing to oxygen or silicon 
as principal binding site for the D X^. An increase in 
the oxygen concentration with increasing iron content is 
confirmed by SIMS data. However, due to the broadening 
of the excitonic transitions the presence of another donor, 
with similar binding and ionization energies cannot be 
ruled out. 

Quenching of the excitonic luminescence upon doping 
with iron is expected, since transition metal impurities 
in III-V semiconductors are well known to act as car- 
rier lifetime killers^ The recovery of the luminescence 
intensity with increasing iron content appears, therefore, 
remarkable and cannot be explained by the increase in 
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FIG. 14: Near-band-edge emission as a function of the Fe 
precursor flow rate. 



shallow donor concentration alone. Apparently, the neu- 
tral donors compete successfully with Fe 2+ in the hole 
capture. The strong quenching of the excitonic emission 
observed at low iron contents seems to be predominantly 
due to the high electron capture rate of Fe 3+ . The change 
of the charge state of iron with the rise of the electron 
concentration in our system is, hence, responsible for the 
recovery of the near band edge emission. 

The recombination of electrons and holes on the 
p e 2+/3+ acce ptor level is accompanied by the character- 
istic intra-center transition of Fe 3+ . Figures [TBI and [T71 
depict the zero phonon lines of the internal, spin forbid- 
den 4 Ti(G)- 6 Ai(S) transition near 1.3 eV, measured at 
20 K under above band-gap excitation in three samples 
with the same iron content (grown at the lowest iron 
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precursor flow rate) but different electron concentration. 



The peaks labelled B, C, D and E at the high energy tail 




as a function of co-doping with Mg. 

Energy [eV] 



FIG. 15: Near band edge photoluminescence in GaN:Fe. The 
positions of D°Xa and the two electron satellite D°XA, n =2 
are indicated by arrows, as well as those of the LO phonon 
replicas for the bound- and free-exciton (FXa). 

It was found that co-doping with Si (Fig. [TTJ) reduces 
the intensity of the Fe 3+ related emission, while co- 
doping with Mg (Fig. [16]) leads to an intensity enhance- 
ment. This behavior is consistent with that of the near 
band edge luminescence. The increased electron concen- 
tration in GaN:Fe,Si as compared to GaN :Fe leads to a re- 
duction of the carrier recombination rate at the Fe 2+ / 3+ 
acceptor level and, thus, the emission intensity, whereas 
reduction of the electron concentration (by codopig with 
Mg) leads to the opposite effect. 
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D. Electron paramagnetic resonance 



FIG. 17: Intensity dependence of the Fe 3+ internal transition 
as a function of co-doping with Si. 



The samples have also been characterized by electron 
paramagnetic resonance. The advantage of this tech- 
nique is that the contributions to magnetization origi- 
nating from different paramagnetic impurities can be de- 
termined separately. It has been found that in all stud- 
ied structures, independent of the iron concentration, an 
EPR signal from isolated, substitutional Fe 3+ ions is ob- 
served. In most samples this signal is superimposed to a 
strong spectrum stemming from some Cr 3+ contamina- 
tion of the sapphire substrate. In other samples (e.g., 
#320) Er 3+ was observed in the substrate instead of 
Cr 3+ . Typical EPR spectra for the magnetic field ori- 
ented along the c-axis are shown in Fig. [T5I 

The spin Hamiltonian parameters for the Fe 3+ ion are 
determined by a least squares fit of the angular depen- 
dence of the EPR line positions for the magnetic field 
B rotated in the (1120) plane - as depicted in Fig. [TH] 



- with the spin Hamiltonian for an S = 5/2 ion in C3,, 
symmetry: 



H = mbB. 9 S - -Bi{Ol + 20^01) - 



B° 2 0° 



BlOl (8) 



Here the first and second term describe the Zeeman- 
and the cubic-crystal-field-interaction respectively, and 
the last two terms represent the second and fourth order 
trigonal crystal field interactions. The spin Hamiltonian 
parameters, especially the value of the second order trig- 
onal parameter I B® , are found to vary slightly with the 
sample thickness, giving indication that the layers are 
strained. 

Typical parameters for a 1 /im thick layer at moderate 
Fe precursor flow rates (50 - 200 seem) are: g\\ — 2.009 ± 
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FIG. 18: EPR spectra - first derivative of the microwave ab- 
sorption with respect to the magnetic field - of Fe 3+ in GaN 
for the magnetic field direction close to the [0001] GaN axis. 
The contribution of Cr 3+ and Er 3+ contamination of the sap- 
phire substrate is easily seen. The lower spectrum was taken 
at a higher microwave power than the upper one (since the 
EPR signal of Er 3+ is not as easily saturated as that of Cr ) 
and a background due to microwave absorption in the cryostat 
is observed. 
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FIG. 19: Angular dependence of the EPR line positions for 
sample #406. The solid lines represent the fit to experimental 
data with the Hamiltonian and parameters given in the text. 



0.002, g x = 2.005 ± 0.002, \B$\ = 237.2 ± 0.6 G, B 4 = 
0.11 ±0.01 G, and \B%\ = 0.8 ±0.1 G. 

Figure [20] shows the concentration of Fe 3+ estimated 
from the EPR signal intensity (circles) compared to the 
total Fe concentration measured by SIMS (triangles) as 
a function of the iron source flow rate for two series of 
samples (full and open symbols). The number of Fe 3+ 



FIG. 20: Concentration of Fe 3+ estimated from the EPR sig- 
nal intensity (circles) as compared to the total Fe concentra- 
tion measured by SIMS (triangles) vs. the iron precursor flow 
rate. Full and open symbols identify two different sample 
series, grown at the same conditions. 



ions in the samples has been determined by compari- 
son with a reference, phosphorus-doped Si sample with a 
known P concentration of 2 x 10 14 cm~ 3 . The measure- 
ments were performed at low microwave powers in order 
to ensure that the intensities of both Fe 3+ and P signals 
were far from the saturation regime. Since the samples 
could not be measured simultaneously, the relative num- 
ber of ions is determined with an accuracy of 20%, the 
uncertainty being due to possible changes in the quality 
factor of the microwave cavity. Moreover, the relative 
concentration contains an additional error related to the 
assumed thickness of the doped layer: the nominal thick- 
ness, obtained from the number of observed thickness os- 
cillations in ellipsometry can differ by 20% from the one 
estimated by SIMS. In the high Fe-doping regime, the ac- 
curacy is even lower, since the surface roughness inhibits 
the observation of the oscillations in ellipsometry as well 
as a reliable estimation of the thickness from SIMS data. 
There is also a systematic error related to the estimation 
of the P content, so that the absolute concentration of 
Fe 3+ determined from the EPR intensity may differ from 
the true value by not more than a factor of two. De- 
spite this uncertainty, a clear trend is visible in Fig. 
with increasing iron precursor flow rate - up to about 
200 seem - the concentration of substitutional Fe 3+ ions 
increases continuously, but at a lower rate than the to- 
tal iron concentration in the sample. Above 200 seem the 
Fe 3+ content is reduced and considerable fluctuations are 
observed from sample to sample. We relate this effect to 
the change of the Fermi level position and, hence, greater 
occupancy of the Fe 2+ charge state. In the high doping 
regime a drastic increase of the oxygen donor concentra- 
tion, by about two orders of magnitude as compared to 
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that of GaN:Fe grown at low Cp2Fe flow rates, is detected 
by SIMS. The oxygen accumulation seems to be corre- 
lated with the enhanced sample roughness: as previously 
underlined, at low flow rates (up to 200 seem) the GaN:Fe 
growth is two-dimensional (and the oxygen concentra- 
tion stays below 10 18 cm~ 3 ), whereas above 200 seem 
the (Ga,Fe)N nucleation is granular with over 1 fim fluc- 
tuations in the grain height at 350 seem. 

The presence of Fe 2+ ions cannot be determined di- 
rectly by EPR since the ground state is a spin singlet 
state and there is no resonance transition observed in the 
X-band. We have detected a metastable 20% increase of 
the EPR signal intensity of Fe 3+ under illumination with 
light of energies above 1.2 eV, confirming the presence of 
Fe in the 2+ charge state, but no estimations of the ac- 
tual concentration can be made from such experiments. 
This increase depends predominantly on the concentra- 
tion of other, not iron-related trap centers in the sample 
as well as on the carrier capture rates. Therefore, in or- 
der to estimate the total concentration of substitutional 
Fe ions in the sample, we have analyzed the EPR line 
width of Fe 3+ . The line shape A(B) of Fe 3+ is found to 
be excellently described by the derivative of the Gaussian 
error function: 



A(B) 



1 



27TCT 



exp 



(B - B )< 
2a 2 



(9) 



where a 2 = (^(B — Bo) ) is the second moment of the 
distribution. At very low concentrations, the lincwidth a 
should be mainly governed by interactions with nuclear 
magnetic moments of the Ga and N isotopes, however, 
with increasing Fe concentration dipole-dipole interac- 
tions between iron ions should lead to line broadening. 
The broadening due to dipole-dipole interaction consists 
of two contributions: one due to the interaction with 
other 3+ ions and one related to the interaction with 
Fe ions in the 2+ charge state. This is true as long as 
the Fermi level is pinned to the Fe 2+ / 3+ acceptor level 
and the oxygen ions have no magnetic moment, which 
seems to be justified, since no EPR transition related to 
O donors is observed. Both contributions depend on the 
i th to ] th ion distance r^ , as r~ 3 . Assuming that the ions 
as well as their charges are uniformly distributed in the 
lattice, one can expect a linear correspondence between 
the dipolar broadening a-ao and the Fe concentration. 
Figure [21] gives the linewidth as a function of iron pre- 
cursor flow rate for the -1/2 — > 1/2 transition of Fe 3+ at 
B || [0001] for the same samples shown already in Fig.lTOl 
As it can be seen, the linewidth increases continuously 
with the Cp2Fe flow rate up to 250 seem and saturates 
above this value. By taking <jq = 12 G for the unbroad- 
ened linewidth we obtain, in the low doping regime, ex- 
actly the same slope for the dipolar line broadening ver- 
sus Cp2Fe flow rate as that for the SIMS concentration 
in Fig. [20] This allows us to set the solubility limit for 
substitutional iron in GaN at our growth conditions (that 
is, the saturation concentration obtained from the EPR 



linewidth) to 1.8 x 10 20 cm" 3 , i.e., 0.4%. It should be 
noted that in contrast to concentration estimations from 
the EPR signal intensity, in the case of EPR linewidth 
there is no error related to the determination of the sam- 
ple thickness, therefore the comparison with SIMS data 
is more straightforward. 
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FIG. 21: EPR linebroadening (left y-axis) as a function of the 
iron precursor flow rate (hollow circles) for <7o = 12 G. The 
SIMS concentration (right y-axis) at low flow rates, denoted 
by full triangles, is shown for comparison. 

The relatively narrow linewidth in the high doping 
regime visible in Fig. [21] indicates that the substitutional 
Fe ions remain well isolated and interact only very weakly 
via dipole-dipole interaction. This is consistent with the 
TEM data, revealing the presence of iron-rich nanocrys- 
tals in which excess iron is accumulated. These nanocrys- 
tals do not contribute to the Fe 3+ EPR. As we already 
underlined, the nanocrystals are located in the upper half 
of the (Ga,Fe)N layer towards the surface and may de- 
plete this part of the layer from isolated Fe ions which 
thus remain in the lower region of the layer only. It is 
worth also noting that the coexistence of Fe 3+ and Fe 2+ 
ions implied by the EPR data as well as by the SQUID 
results, suggests that the Fermi level is pinned well below 
the conduction band minimum by the Fe 3+ /Fe 2+ accep- 
tor states. Such pinning is in contradiction with the rel- 
atively high value of electron concentration and mobility 
we observe in this high Cp2Fe flux range, as shown in 
Figs. [TO] and [TOJ giving clear indication that the Fermi 
level is located either in the conduction band or in the 
donor impurity band. The simultaneous observation of 
Fe 3+ ions (as evidenced by EPR) and free electrons (as 
shown by the Hall data) in the high Fe concentration 
regime is one of the surprising findings of the present 
work. In view of our TEM data, it is natural to suggest 
that the conduction proceeds via the upper region of the 
layer depleted from isolated Fe acceptors, while Fe 3+ ions 
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exist mainly in the part of the layer closer to the interface 
with the GaN buffer. 



E. Magnetic properties of (Ga,Fe)N 

1. Paramagnetic contributions 

Figure [22] gives the temperature dependence of the 
magnetic susceptibility x(T) obtained from SQUID mag- 
netization measurements at a constant magnetic field 
of 0.1 T for four GaN:Fe samples with a low Fe con- 
tent. As it can be seen, two components of x(T) may 
be distinguished: one is a Curie-type of paramagnetism 
X(T) ~ l/T, which we assign to Fe 3+ (d 5 ) ions. The con- 
centration of Fe 3+ ions obtained under this assumption is 
shown in Table I for the samples under investigation. The 
second component is a temperature independent contri- 
bution dominating at high temperatures. However, as 
shown in Fig.l23( the magnetization at high temperature 
in addition to the paramagnetic part linear in the mag- 
netic field, contains a ferromagnetic component. Proper- 
ties and origin of the high temperature ferromagnetism 
in (Ga,Fe)N will be discussed later. 
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FIG. 22: Magnetic susceptibility as a function of tempera- 
ture, showing the temperature independent and the temper- 
ature dependent components assigned to the Van Vleck and 
Curie paramagnetism, respectively. The ferromagnetic and 
substrate contributions (temperature independent and linear 
in the field) have been subtracted. 

We attribute the term linear in the magnetic field to 
Van Vleck paramagnetism, most likely originating from 
the Fe ions in the 2+(d 6 ) charge state. The presence 
of Fe ions in such a charge state seems to correlate with 
the Hall effect measurements which, according to Fig. [TTJ 
indicate that the Fermi level for low Fe flow rates may be 
pinned to the Fe 3+ / 2+ level. The coexistence of Fe in the 
3+ and 2+ charge states is also confirmed by photo-EPR 
results and the theoretical explanation of the Van Vleck 
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FIG. 23: Magnetization from SQUID measurements at 200 K 
for a series of (Ga,Fe)N samples with different Fe content 
(solid symbols). The substrate contribution (temperature in- 
dependent and linear in the field) has been subtracted in this 
and in the figures that follow. Data obtained by numerical 
reflection, if only a half of the hysteresis cycle has been mea- 
sured, are shown by open symbols. 



contribution is currently being formulated^ 



2. Ferromagnetic properties 

Ferromagnetic signatures, such as magnetic hysteresis, 
remanent magnetization (M rcm ), and coercivity (H c ) are 
detectable in all our Fe-doped films but with the increase 
of the total Fe concentration above the solubility limit of 
0.4% the ferromagnetic component exceeds the paramag- 
netic contributions previously discussed. This is depicted 
in Fig. [24] which shows the magnetization cycles at 5 K 
for samples with various Fe contents. A rather clear ferro- 
magnetic component is seen, with its strength increasing 
upon the increase of the Fe content. For the layers with 
the highest Fe concentration, we recover about 25% of 
the Fe ions contributing to the saturation value of mag- 
netization, if S = 2 and g = 2 are assumed. 

According to the findings summarized in Fig. [25] the 
ferromagnetic signal does not diminish remarkably when 
the temperature is increased to 200 K, and it persists up 
to above room temperature, as shown in the inset. The 
set of magnetization curves obtained in the full avail- 
able temperature and field range allows to establish the 
temperature variation of spontaneous magnetization Ms, 
determined from the Arrot plot, as depicted in Fig. [26] 
If the ideal Brillouin Ms(T) dependence is assumed, the 
data displayed in the inset point to an apparent Curie 
temperature Tq of ~ 500 K. The ferromagnetism in ques- 
tion reveals a well defined anisotropy. Figure [27] shows 
that the saturation value of magnetization is greater for 
the magnetic field perpendicular to the wurzite c-axis. 
Interestingly, a similar anisotropy was observed in the 
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FIG. 24: Magnetization at 5 K for a series of (Ga,Fe)N sam- 
ples with various Fe contents (solid symbols). Data obtained 
by numerical reflection, if only a half of the hysteresis cycle 
has been measured, are shown by open symbols. 
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FIG. 25: Magnetization at 200 K for a series of (Ga,Fe)N 
samples with various Fe content (solid symbols). Inset: above 
room temperature M(H) for sample #408. Data obtained by 
numerical reflection, if only a half of the hysteresis cycle has 
been measured, are shown by open symbols. 



case of (Zn,Co)0^ 

We want to underline here that we do not take all 
the above ferromagnetic signatures as a solid proof of a 
uniform ferromagnetic state of (Ga,Fe)N. In the mate- 
rial under consideration here, the average concentration 
of magnetic ions is far below the percolation limit for 
the nearest-neighbor coupling, and at the same time the 
free-carrier density is too low to mediate an efficient long- 
range exchange interaction. Hence, we take the presence 
of high temperature ferromagnetism as an indication of 
the non-uniform distribution of the Fe ions over the GaN 



FIG. 26: Determination of the spontaneous magnetization 
at various temperatures obtained by plotting the square of 
the magnetization vs. the ratio of the magnetic field to the 
magnetization (Arrot plot, sample #408). The inset shows 
the established spontaneous magnetization as a function of 
temperature, whose extrapolation leads to the apparent Curie 
temperature of 500 K for the same sample. 



lattice already demonstrated via TEM analysis of the 
structures. 

In order to shed some light on the origin of the ferro- 
magnetic behavior presented in the last figures, we dis- 
cuss at first the results from zero-field cooled (ZFC) and 
field cooled (FC) magnetization measurements at a low 
magnetic field of 50 Oe, given in Fig. [55] The markedly 
different behavior, namely the increasing FC-ZFC differ- 
ence on lowering temperature, is indicative of the fact 
that the sample is in a magnetically frozen state hav- 
ing its response determined by thermally activated pro- 
cesses across some energy barriers. We further note that 
the ZFC and the FC magnetizations are characterized by 
temperature gradients of opposite sign without a maxi- 
mum appearing in the ZFC curve. Such a behavior can 
be expected for an ensemble of magnetically anisotropic 
ferromagnetic particles thermally cycled below its block- 
ing temperature (Tb). At T < Tb the moments of some 
of the particles are blocked and unable to surmount the 
magnetic energy barriers at the time scale of the exper- 
iment. The energy barriers in the vicinity of H = 0, 
cannot be overcome without an increase in the thermal 
energy (or an external magnetic field). For a system of 
non-interacting, uniaxial particles the energy barrier for 
the magnetization reversal is AE — KV, where V is 
the particle volume and K is the anisotropy energy den- 
sity. Therefore, the actual magnetic behavior depends 
on the measuring time (t m ) of the specific experimental 
technique with respect to the time t required to over- 
come the energy barrier (for DC magnetic studies a t m 
of 100 s is usually adopted). Accordingly, for a system of 
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FIG. 27: Magnetization measurements for sample #408 (solid 
symbols), indicating that the saturation value of magnetiza- 
tion is greater for the magnetic field Hu along the surface of 
the sample (perpendicular to the c-axis). Data obtained by 
numerical reflection, if only a half of the hysteresis cycle has 
been measured, are shown by open symbols. 



This is a rather high value and it indicates that the fer- 
romagnetic particles in (Ga,Fc)N contributing to the fer- 
romagnetic response must be characterized by substan- 
tial values of the KV product. The existence of fine en- 
tities in the studied samples is corroborated by our TEM 
analysis, revealing the presence of wurzite nanocrystals 
embedded in the host GaN. The substantial amount of 
these nanocrystals, as measured by the value of Ms, is 
clearly triggered by higher flow rates (above ~150 seem) 
- see Fig. [221- and undoubtedly correlates with the nom- 
inal Fe concentration in the layers. Furthermore, as al- 
ready pointed out, EDS results demonstrate that the Fe 
concentration is largely enhanced within the nanocrys- 
tals. 
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FIG. 28: Magnetization of sample #466 measured at H appL = 
50 Oe after cooling the sample in absence of the external field 
(ZFC) and during cooling in the same H appL . Qualitatively 
the same results are obtained at H appL = 200 Oe. 



identical non-interacting uniaxial particles, the following 
simple condition for the evaluation of Tb is employed: 



T B ~ KV/25k B . 



(10) 



Therefore, T > Tb defines the onset of superparamag- 
netic behavior for this system, whereas below Tb an ir- 
reversible behavior is expected. Within this line of inter- 
pretation, our experiment (Figs. |2"51|2"8)) points then to a 
large, at least above room temperature Tb for (Ga,Fe)N 
layers. 
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FIG. 29: Saturation magnetization of the ferromagnetic com- 
ponent extracted from M(H) curves at 200 K for the #466 - 
#473 layer series versus the Fe precursor flow rate. 



All these findings imply that in the case of Fe con- 
tents above the solubility limit (Ga,Fe)N undergoes a 
spinodal decomposition into regions with small and large 
concentrations of magnetic ions, so that most of Fe ions 
resides in the magnetic nanocrystals. Now, taking the 
volume of these Fe-rich nanocrystals as V — 1000 nm 3 
on average, we end up with a requirement for K to be 
~ 20 x 10 5 erg/cm 3 to reproduce a T B of about 500 K, 
as the magnetization studies indicate. This value is four 
times higher than that of bulk cubic iron and it is a half 
of that of hexagonal cobalt. However, the estimated large 
value for K should be treated rather as an upper limit, 
since Tb can get enlarged above the value predicted by 
Eq. [TO] if, as a result of the magnetic interaction between 
particles, an ordering of the magnetic moments takes 
place. Moreover, we note that the nanocrystals are not 
monodispersed and that the distribution in their volumes 
gives rise to a distribution in the blocking temperatures. 
Thus, Eq.fTOl gives the mean Tb for the system, so that ir- 
reversible effects occur also above Tb and, therefore, the 
requirements for the magnitude of K are considerably 
relaxed. 
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V. SUMMARY AND OUTLOOK 

Our results on the GaN:Fe and (Ga,Fe)N systems have 
reemphasized the importance of two related and generic 
aspects of magnetically doped wide band-gap semicon- 
ductors. First, most of the substitutional magnetic im- 
purities give rise to the presence of deep levels originat- 
ing from open d or / shells of the impurity atoms. The 
presence of such levels strongly affects the electrical and 
the optical properties of these materials systems. At the 
same time, co-doping with shallow impurities, by chang- 
ing the charge and, thus, the spin state of the magnetic 
ions, can serve to alter the magnetic response. Second, 
though depending on the growth conditions, the TM sol- 
ubility limit is rather low and typically - except for Mn - 
it is hard to introduce more than 1% of magnetic im- 
purities into randomly distributed substitutional sites. 
Due to the low density and in absence of band carriers, 
the impurity spins can be regarded as decoupled. Ac- 
cordingly, pertinent properties of the system can be ad- 
equately described by the single-impurity theory^ The 
series of transport, optical, EPR, and magnetic exper- 
iments that we have reported here for MOCVD grown 
(Ga,Fe)N and for (Ga,Fe)N co-doped with either Si or 
Mg has substantiated the notion that the material char- 
acteristics are determined by the Fe 3+ /Fe 2+ acceptor- 
like deep level. Thus, the relative concentrations of Fe 3+ 
and Fe 2+ ions and, hence, the relative importance of the 
Curie and Van Vleck paramagnetism is defined by the 
density of shallow acceptors and donors. This means, in 
particular, that if the acceptor concentration were high 
enough to give rise to the presence of a sufficiently in- 
creased density of weakly localized or delocalized valence 
band holes, the carrier-mediated coupling between spins 
localized on Fe 3+ ions would operated Our experiments 
have indicated also that the solubility limit is about 0.4% 
in the case of Fe in GaN under our growth conditions. 

A number of scenarios can be realized when increasing 
the Fe concentration beyond the solubility limit. One 
would be the appearance of Fe in interstitial positions, 
where it will act as a double donor. We have indeed 
found a strong increase in the electron concentration in 
this regime. EXAFS or Rutherford backscattering exper- 
iments might allow to detect the presence of interstitial 
Fe. A further scenario is represented by the crystallo- 
graphic phase separation in the form of precipitates of 
Fe or Fe-compounds. Finally, an appealing possibility is 
the existence of spinodal decomposition into regions with 
respectively high and low concentrations of the magnetic 
constituent. Our TEM results have demonstrated the 
presence of coherent Fe-rich nanocrystals embedded and 
stabilized by the GaN host. Furthermore, the TEM imag- 
ing shows that these nanocrystals are not distributed uni- 
formly over the layer but float towards the surface during 
the growth process. The resulting inhomogeneity in the 
layer structure along the growth direction may account 
for the increased surface roughness at high doping level, 
as observed by AFM. This inhomogeneity may also give 



reason for the simultaneous presence of deep Fe acceptors 
and conducting electrons demonstrated by our transport, 
EPR, and SQUID experiments. We argue that the Fe- 
rich nanocrystals account for the high-temperature fer- 
romagnetic properties of the (Ga,Fe)N films grown by 
MOCVD and MBE as well as of GaN implanted with 
Fe. We note that the GaN host may stabilize nanocrys- 
tals of a composition and structure non-existing in the 
case of a free standing Ga-Fe-N material systems. Actu- 
ally, such a situation has lately been found in the case 
of MBE grown (Ge,Mn), where spinodal decomposition 
into Ge and novel Mn-rich (Ge,Mn) nanowires has been 
revealed4£ Furthermore, it has recently been suggested 
by one of us^ that the aggregation of magnetic nanocrys- 
tals can be modified by altering the charge state of mag- 
netic impurities during the growth or by thermal treat- 
ment during or upon deposition. This can obviously be 
accomplished by co-doping with shallow donors and ac- 
ceptors. Experiments in this direction are on the way. 
While our paper is concerned with one material system, 
namely GaN:Fe - (Ga,Fe)N, there is ground to conjecture 
that spinodal decomposition and inhomogeneity along 
the growth direction revealed by our work, as well their 
influence on transport, optical, and magnetic character- 
istics, are generic properties of diluted magnetic semicon- 
ductors and diluted magnetic oxides showing high appar- 
ent Curie temperatures. Within this model, the sponta- 
neous magnetization originates from magnetic moments 
of ferromagnetic or ferrimagnetic nanoclusters or from 
uncompensated spins at the surface of antiferromagnetic 
nanoparticles. We are confident that the identification 
of the actual chemical nature of nanocrystals in partic- 
ular matrices and the determination of the mechanisms 
accounting for the elevated blocking temperatures nec- 
essary to make the ferromagnetic-like characteristics to 
survive up to high temperatures, will attract considerable 
attention in the years to come. 
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Appendix: Experimental difficulties and precautions 
in SQUID measurements 

General remarks - the (Ga,Fe)N layers are grown on 
330 /im thick sapphire substrates, so they constitute only 
a fraction of the whole volume of the investigated mate- 
rial, and their magnetic moment is very-small to small 
when compared to the diamagnetic signal of the sub- 
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strate. Since it is unfeasible to remove the sapphire 
substrate, the magnetic moment of the (Ga,Fe)N system 
has to be established by post-measurement subtraction of 
the known magnetic moment of a Fe-free, but otherwise 
identical GaN/sapphire reference sample. This in general 
does not pose any experimental difficulties if the signal 
to be investigated is at least comparable to the substrate 
signal. In the case of our epitaxial layers of magnetically 
diluted materials, the weak magnetic response, as shown 
later, requires typically more than 90% compensation, 
approaching 99% for the lowest Fe concentration. Fur- 
thermore, the overall measured signals are rather weak, 
so there is little room left for a significant improvement 
of the absolute accuracy of the experimental data af- 
ter the compensation, as we are approaching the base 
noise level of the magnetometer set up. Results of typ- 
ical measurements are depicted in Fig. 1301 exemplifying 
the necessity for high quality, distortion free measure- 
ments. Additionally, any experimental flaws, short-cut 
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FIG. 30: Typical examples of magnetic moment vs. magnetic 
field m(H) recorded for the layers presented in this paper. 
Left panels: raw data. Right panels: after subtracting the 
major linear component. The high slope lines on the left pan- 
els indicate the uncompensated m(H) data. Data obtained 
by numerical reflection, if only a half of the hysteresis cycle 
has been measured, are shown by open symbols. 



procedures or negligences in the sample handling have se- 
vere adverse effects on the reliability of the final results, 
as they get strongly magnified by the high degree of com- 
pensation analysis required to determine the signal from 
these layers. To preserve the credibility of the measure- 
ment, tight constrains on the whole measurement process 
are imposed and can be generalized as: contamination- 
free handling of the samples, magnetically-clean mount- 
ing in the SQUID magnetometer, and credibility in the 
determination of the magnetic moment by the magne- 
tometer software. 



a. Handling of the samples We perform all the man- 
ual operations, including the sample mounting on the 
magnetometer sample holder, in a laminar flow cabi- 
net equipped with an absolute filter. Also, we avoid to 
bring the samples into contact with any ferrous mate- 
rials, including metal tweezers^ 5 - since (stainless) steel 
tweezers are a well recognized source of contamination 66 
of, e.g., silicon wafers and their magnetic signature can 
be comparable with the magnetic response expected in 
our studies. 65 However, contrary to Si, where the con- 
tamination can occur either by contact chemical reac- 
tion or by mechanical scratching, we believe that in 
our case mostly a mechanical grinding of the softer 
metal by the much harder sapphire/GaN is the major 
source of magnetic contamination. Therefore, we employ 
only titanium-/plastic-based tweezers. Another poten- 
tial source of specimen contamination is related to the 
dicing of the sapphire wafer, as it requires a diamond 
saw. For the preparation of our samples, we employ only 
plastic-based cutting wheels on the wafer coated with 
photoresist for additional protection. After dicing, the 
samples are rinsed with organic solvents in an ultrasonic 
cleaner. We verified that none of these substances carries 
any detectable amounts of ferrous metals. 

b. Sample holder and sample mounting In the 
present study we utilize a home-made SQUID magne- 
tometer built into a standard cryostat operating from 5 
to 325 K in continuous flow mode. Although equipped 
with a 1 T superconducting coil, the set-up has been op- 
timized for very weak field operation and the noise level, 
increasing with the field, limits our experimental abili- 
ties to ~ 0.5 T. We use a commercially available SQUID 
sensor and the matched electronics, connected to a com- 
puter based acquisition system running under an home 
developed code. The sensing end of the set-up, like most 
modern similar systems, utilizes the controlled tripping 
of the sample (with the sample holder) through sensing 
coils arranged in the 2nd order gradiometer configura- 
tion. Ideally, only the sample to be measured produces 
the signal. Its position-related voltages generated by the 
SQUID electronics during this movement form the base 
for the numerical procedures necessary for extraction of 
the magnetic moment. In general, a gradiometer config- 
uration allows the use of any uniform extended object - 
rod-like shaped - as a sample holder, since if its length 
is much larger than the extent of the pick-up coils, the 
sample holder produces no output according to symme- 
try considerations. In reality, the limited length of the 
sample chamber does not allow sufficiently long sample 
holders, resulting in a non-linear position-dependent volt- 
age disturbing the signal from the sample, particularly 
in the most sensitive ranges. This parasitic contribution 
can be largely reduced by employing magnetically weak 
substances, preferably diamagnets, shaped in the form of 
small tubes or rods. Nothing, however, will compensate 
for the shape variations and for the distribution of in- 
homogeneous magnetic impurity inside the rod material. 
The former result from all sort of mechanical imperfec- 
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tions (e.g. variable cross section area, scratches, etc.), 
the latter may be an inherent property of the chosen ma- 
terial. Both generate position-dependent voltages that 
modulate or distort the signal from the sample in an un- 
known and a-priori unpredictable T- and H- dependent 
way, corrupting, in this way, the results of the measure- 
ments in an unrecoverable way. Therefore, the testing of 
the sample holder and the selection against the parasitic 
magnetic signal turns out to be of a great importance for 
this study. 

Currently, the most popular choice for magnetometer 
sample holders are drinking straws, not only because of 
their low magnetic signature, but also due to the flexi- 
bility allowed in a glue-free mounting of samples with di- 
mensions comparable with the diameter of the straw. For 
the purposes of our studies, we have undertaken an ac- 
curate testing of the available materials, focusing in par- 
ticular on those directly provided by Quantum Design. 67 
Some 30 straws were scanned along their length at 3000, 
1000 Oe, and at remanence, usually at low T (10 K), 
and occasionally at 200 K. All tested straws have been 
thoroughly cleaned in organic solvent prior to the tests. 
The collected position dependent signals are processed 
according to the same numerical steps employed to estab- 
lish the magnetic moment in the magnetometer, so that 
the spatial distribution of the magnetic moment along 
the straws could be established. We have found that al- 
though on average the straws have a small magnetic sig- 
nature suitable for general magnetic investigations (ex- 
hibiting moment variations below 2xl0~ 7 emu on a one 
inch span near their central part), only two went below 
~ 5xl0~ 8 emu, reaching in this way the upper limit re- 
quired for our studies. Instead of using potentially prob- 
lematic straws, we opted for ~ 20 cm long and 1 or 2 mm 
wide Si strips cut from 8" commercial wafers^ After a 
thorough testing, we have found that all the strips con- 
form to our stringent requirements concerning the mag- 
netic uniformity. Additionally, the Si strips ensure an 
excellent rigidity and allow for an easy sample mounting 
in various in-plane orientations. The major drawback of 
switching to Si strips is the necessity of using a glue to fix 
the sample onto the strip. After testing many substances 
for a possible contribution to the measurement as well as 
for mechanical reliability and low temperature behavior 
we settled for the heavily diluted GE varnish, which best 



fulfills all requirements. We found that only the recently 
available orange version^ as opposed to the dark red 
one, shows no magnetic signature down to 5 K and up 
to 3000 Oe even when examining ~10 times the amount 
normally needed for mounting. It must be noted, that the 
usage of narrow strips may lead to larger sample off cen- 
tering than in the case of wider straws. In our case how- 
ever, owing to the relatively large diameter of our pick-up 
coils we determined this effect to be below 0.5%^ i.e. 
within the typical run-to-run non-reproducibility of our 
system. 

c. Reduction of the data Typically the investigated 
magnetic moment is established by a numerical fitting 
(least square regression) of an ideal, point dipole re- 
sponse function of the gradiometer pickup coils into the 
set of voltages collected during the tripping of the sam- 
ple through the coils. In the case of macroscopically large 
samples of rectangular shape the form of the real response 
function starts to differ noticeably if the specimens oc- 
cupy more than a few percent of the volume of the pickup 
coilsj 7 ^ As a consequence, even the best possible fitting 
routine will not establish the accurate sample moment, 
unless the response function gets adequately modified. 
It is relatively easy to extend the point dipole response 
function to account for uniformly magnetized body of a 
given lengt h 70 i 71 i 72 and this length parameterized func- 
tion is routinely used in our system to extract a length- 
independent magnetic moment. This however seems to 
be the only one feasible to compensate sample size ef- 
fects. The account of other detrimental effects like non 
uniform sample extension in the plane perpendicular to 
the magnetometer axis or displacement of the whole sam- 
ple aside from the axis is too complex to get generalized 
in any sort of universally parameterized function. So, in 
order to minimize their adverse influence on the experi- 
mental results, we confined ourselves to one sample size 
for magnetic measurement and additional care is taken to 
mount the samples on Si sticks in a reproducible way. For 
completeness, we note that from the results of numerical 
modelling 70 ' 71 a maximum 1% of moment distortion is 
expected for the samples used in the experiment. There- 
fore, we can safely conclude that although potentially 
large, the sample-size-related effects should not system- 
atically alter the magnetometer results before and after 
the compensation. 
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